arXiv:hep-ph/9706328v3 8 Oct 1997 


Production process dependence of neutrino flavor conversion 


Subhendra Mohanty 

Theory Group, Physical Research Laboratory, 
Ahmedabad - 380 009, India 


Abstract 


We perform a covariant wave-packet analysis of neutrino oscillations tak¬ 
ing into account the lifetime of the neutrino production process . We find 
that flavor oscillations in space are washed out when the neutrinos are pro¬ 
duced from long lived resonances - and what may be observed in appear¬ 
ance/disappearance experiments is a uniform conversion probability indepen¬ 
dent of distance. The lifetime of the resonance which produces the neutrinos 
acts as the the effective baseline of the experiment. For this reason the LSND 
experiment where neutrinos are produced from muon decay has two orders 
of magnitude more sensitivity to neutrino mass square difference than other 
experiments where the neutrinos are produced from pion or kaon decays. 
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There are two different formulas which describe flavor oscillations in spacetime; the one 
applied to neutrino oscillations [1] is derived in the high energy regime and the other applied 
to kaon oscillations [2] is valid in the non-relativistic limit. From the derivations of these 
formulas it is not clear what formula is applicable at intermediate energies (for example for 
kaon oscillations from stopped protons where {m ^/~ 25%) . A covariant derivation of 
a flavor oscillation formula which would be valid at both high and low energies would be of 
interest from the conceptual as well as the experimental point of view. 

Kayser and Stodolsky [3] have advocated a covariant generalization of the non-relativistic 
phase factor exp{—imt) by replacing the absolute time t by the Lorenz invariant proper 
time s = {t^ — In the lab frame the covariant expression for the phase factor is 

exp — irriiSi = exp {—imft/Ei). The phase difference between two mass eigenstates is 
/E) ~ (Am'^/E) which is twice the the phase difference of the standard formula [1]. 
This result has prompted the claim [4] that kaon oscillations in <F factory will have oscillation 
length which is half of what is given by the standard formula. This claim has been refuted 
by [5] who take the view that the interference phase difference should not be evaluated at 
different space-time points but should be evaluated at the average spacetime interval and 
taking the phase difference to be (mi — m 2 )(<si -|-S 2)/2 instead of (miSi — 171282 ) they recover 
the standard formula. Other methods of showing that an extra factor of two does not 
appear in the kaon oscillations formula have been discussed in [6]. A covariant derivation 
of the neutrino oscillation formula has been given by Grimus and Stockinger [7] who treat 
the entire process of neutrino production ,propagation and detection as a single Feynman 
diagram. They show that on taking the large distance limit of the neutrino propagator the 
scattering cross section shows a space-time oscillatory behaviour and the oscillation length is 
identical to that given by the standard formula [1]. G-S assume the initial states to be plane 
waves therefore the concept of a coherence length [8] does not emerge in their formulation. 

In this paper we derive the oscillation amplitude by evaluating the Feynman propagator 
of the neutrinos in the large time-like asymptotic limit. The asymptotic propagator of 
a position eigenstate has the form K{x,t;mi) ~ {mi/27risi)^^'^ exp—irriiSi where Si is the 
spacetime interval propagated by the m* mass eigenstate. This is a held theoretic derivation 
of the Kayser-Stodolsky [3] phase factor. If the initial wave-functions of the propagating 
particles were strictly delta functions in spacetime then no interference between different 
mass eigenstates can take place. We therefore generalize the delta function propagators to 
propagators of Gaussian wave-packets. The interference term as a function of distance is 
obtained by taking the time-overlap of different mass eigenstate propagators . The expression 
for havor conversion probability (for say two havors with mixing angle 6) as a function of 
distance X turns out to be, 

1 /\ iry^ 2 

P{a ^ 7, X) = j W29(l - cos{—X) (1) 

Therefore extra factor of two which came from naively subtracting the phases of plane wave 
propagators [4] goes away in the wave-packet averaging and the standard expression for the 
oscillations length is recovered. This is in agreement with [5]- [6]. The covariant wave- 
packet treatment introduces a new contribution to the exponential factor A when particles 
are produced from long lived resonances ( for example for neutrinos from muon decay as 
opposed to neutrinos from Z decay). At distances smaller than the coherence length , the 
exponential suppression factor is, 
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A 


, Am^ 
'2V2B 
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where r is the lifetime of the resonance which produces the particle which undergo flavor 
oscillations. If the uncertainty in position of the initial particle vr is larger than the detector 
distance X , the exponential term washes out the oscillations in (1). In neutrino experiments 
where the source is pions, kaons, muons or nuclear fission, the spatial oscillations of the 
conversion probability cannot be observed. What can be observed in these experiments 
is a constant (distance independent) conversion probability. The conversion probability is 
sensitive to values of Am?' ~ {2^/2E/t) , which means that the experimental bound on 
Am^ is lower with longer -lifetime sources. For this reason the LSND experiment [9] which 
uses neutrinos from muon decay (r^ = 2.19 x 10“®s) is sensitive to two orders of magnitude 
lower neutrino mass square difference compared to other accelerator experiments like BNL- 
E776 [10], Karmen [11] and CCFR [12] which use neutrinos from pion and kaon decays 
(r ~ 10“®sec). We £t the experimental data from LSND along with BNL-E776, Karmen, 
CCFR and Bugey [13] experiments with the conversion probability formula (1) and plot the 
allowed range for the mass difference and mixing angle. The lower bound on allowed Am^ 
is two orders lower than what is obtained by fitting the LSND muon decay data with the 
standard oscillation formula (equation (1) with A set to zero). 

Field theoretic derivation of phase factor: The phenomenon of flavor oscillations takes 
place because particles are produced and detected as weak interaction eigenstates (the neu¬ 
trino states z/g, Vfj,, Vt or the Kaon states 77°, etc) but the propagators are diagonal in the 
mass eigenstates (the neutrino mass eigenstates z/j,i = 1 — 3 or the Kaon mass eigenstates 
Kl, Ks). The probability amplitude for oscillations of a gauge eigenstate (|q; >) to another 
(1/5 >) is a linear superposition of the propagation amplitude of the mass eigenstates (|i >). 


A{a —> /5;t) = ^ < j3\i> < i\e ^^^\i >< i\a > 

i 


( 3 ) 


In relativistic field theory the propagation amplitude of the mass eigenstates < > 

can be identified with the Feynman propagator Spixf — Xi^m? =< T Ui{xi)t)i{xf) >. The 
Feynman propagator in position space is , 


K{x, m 



d^P ^ + 

( 27 r)^p^ — mf —ie 


( 4 ) 


This integration can be done by expressing the denominator as an exponential. 


—i 


p^ — — ie 


roo 

/ da exp {ia{p‘^ — — ie)} 

Jo 


( 5 ) 


and integrating over the resulting Gaussian in p. The remaining integral over a, 

,.(t^-|xp)- 


1 -2 r- 2 -R - 

/ da a exp {lam, -1 

Jo ^ ^ t 4q; 


( 6 ) 


can be performed by substituting s = {t^ 


^)^/^, = ims and rj = 2(am/s) and making 


use of the integral formula [14] for the Bessel function 
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( 7 ) 


1 7 ^ 

- dr] exp - -(77 + -) = K^{z). 

2 Jo 2 p 

The resulting expression for the Feynman propagator (4) is 

7 7Tl • 

K{x,mi) = {^){i P + rrii) (—) Ki{imiSi) ( 8 ) 

47r^ Si 

where s* = is the invariant spacetime interval propagated by the z/j mass eigen¬ 

state. If this interval is large (s >> m~^) and time-like ( t > |x| ) then we can use the 
asymptotic expansion of the Bessel function 

/ 2 

KAims) ~ (-) e (9) 

\mmsj 

to obtain from (8) the expression for the propagation amplitude at large time-like separation 

K{x,t;mi)=( - ^ \ exp{-im^{tf - tf - (xf - Xi2)} (10) 

\2m^{tf-ty - (xf -Xi)2)y 

The phase factor of the amplitude (10) is Lorentz-invariant and can be written in terms of 
the neutrino energy E and the time of flight t as measured from the lab frame as —irriiSi = 
—irrii (1 — t = —iirrif /Ei) t. 

An extra factor of two appears on subtracting the phases at different spacetime points. 
A flavor eigenstates neutrino or kaon is observed at a single spacetime point and one should 
therefore compute the phase difference at the overlap of the two mass eigenstate wave- 
packets. This averaging over spacetime is done formally by considering the propagators of 
Gaussian wave-functions as opposed to plane waves , and the standard expression for the 
oscillation length is recovered. 

Conservation laws for long distance propagators: A flavor eigenstate propagates as a 
linear combination of different mass eigenstates which are on-shell. The phase difference 
is obtained by some authors by assuming that the different mass eigenstates have common 
energy, and by some by assuming that they have a common momentum [15]. In this section 
we show that the linear combination of different mass eigenstates have neither the same 
energy nor the same momentum. We show that particles propagating over large distances 
{X » P/mf) are constrained by the conservation laws at the vertex to be on shell. The 
on-shell condition is all that is needed to fix the phase difference and the oscillation length. 

Consider a diagram with a propagator between vertices at spacetime points (xi, 0 : 2 ) with 
a number of external legs at the vertices. The amplitude is proportional to 

J d^xi d‘^X2 exp{—i'^qi-xi) K{xi,X2) exp{i E?/-^2) (11) 

* / 

Where g* are the incoming four momenta from the external legs at xi and g/ are the incoming 
four mommenta from the external legs at X2- Substituting the propagator 

/ g-ip-{xi-X 2 ) 

d^P— -y—— (12) 

p^ — + le 
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in (11) and integrating over xi and X 2 gives the usual energy mommentum conservation 
<ii — p) and g/ — p) at the vertices. When the proper distance between the two 

vertices s is larger than m~^, then the propagator used in (11) is the asymptotic form (10), 
which in mommentum space can be written as [16], 

K{xi,X2;s » m-^) = 

= / /-ot ^ 2 {-i(VP2 + m2)(t/ - U) + iP ■ (xf - Xi)} (13) 

J V P + 

Substituting (13) in (11) and integrating over xi we obtain 

J d!^xi exp {-f (XI - VP^ + w? )ti + i(X Qi - P) ' xi} 


= VP2 + m2) <53(Xqi-P) (14) 

i i 

This implies that at a vertex energy and mommentum are conserved and in addition those 
particles which propagate without freely over distances larger than (|P|/m^) obey the mass 
shell constraint , E = 

Wave-packet analysis: The propagator K{xf,Xi) given in (10 ) is the amplitude for 
the propagation of a particle localized at Xi (a delta function initial wave function) to be 
detected at Xf. The propagator of some general initial wave-packet 'hi„(a; — Xi) is obtained 
by using the expressions for the delta function propagator given in (10) and the superposition 
principle, 

K{xf-Xi)= J d‘^x K{xf,x) ^in{x - Xi) (15) 


If the initial wave function is a Gaussian , 

'^inix-Xi) = N exp{-iEa{t-ti)+iP^-exp - 2 '^ } (16) 

where and at are the uncertainties is the initial position and time of production of 
the particle respectively and N is the normalization constant. In the earlier wave-packet 
analyses of oscillation problem [8], the uncertainty in the time of production was neglected. 
In a covariant treatment both should be included. In the next section we show that the 
time uncertainty gives rise to a novel phenomenon of conversion without oscillations in 
many experimentally relevant situations. Substituting the expression for K{xf,x) given in 
(10) and the Gaussian initial wave-packet (16) in (15) and evaluating the integral by the 
stationary phase method [17] we obtain the expression for the propagation amplitude of a 
Gaussian wave-packet , 


K{X,T;ma) = 


■47r|X| 




n 




^-1/4 


.m 


exp{-t-^T +*Pa-(X-VaT)} 


X 


exp{ 


(X-Va T)\ 


(17) 
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where X = Xf — Xi and T = tf — ti are the space and time intervals propagated by the 
center of the Gaussian wave-packet and Va = Pa /Ea . The normalization constant has been 
hxed such that / dT K(X.,T;ma) = 1. The factor of (l/ 47 r|Xp)^/^ in 

(17) shows that the particle flux decreases inversely with the sqaure of distance. Neutrino 
disappearance experiments look for evidence of depletion of a certain neutrino species over 
and above the expected inverse square decrease in flux. In the following we will not display 
this factor (l/ 47 r|Xp) in the probability expressions. The amplitude for the oscillation of one 
type of neutrino flavor to another is obtained by substituting the propagation amplitudes 
for mass eigenstates (17) in , 

3 

A{a^p-X,T) = Y.Uo.,K{ma]X,T) (18) 

a 

where a, (3 denote the flavor eigenstates (z/g, or K^,K^) and the summation index a 

denotes a mass eigenstates (z/i, z/ 2 , z /3 or Kl, Ks) ■ Uaa =< o:\a > and =< a\P > are the 
elements of the mixing matrix which relate the flavor eigenstates with the mass eigenstates. 
The probabilty of flavor oscilation as a function of space-time is the modulus squared of the 
amplitude (27) 


Pia^/3;X,T) = \J2U^aKima;X,T) U, 


* 12 
af3 1 


= J2\UfSa\" \Uaa\" + T. ^fSa K,U;, K{X,T;ma) K^X, T; TTlfe) (19) 

•T a^b 

In interference experiments the time of flight of the particle is not measured and only 
the distance between the source and the detector is accurately known [15] . The probability 
for the flavor conversion as a function of distance is given by the time integral of (19) , 


P{a^(3-X)= / dT \ K{ma,X,T) 


= Y.\Upa? \Uc.a? + [ dTK{X,T-ma) K\X,T-m,) {20) 

O' a^b 


The interference term given by the time-overlap of the propagation amplitudes (wave- 
functions) of different mass eigenstates can be evaluated for the Gaussian propagator (17) 
and is given by, 


Re [dT k{X,T-ma) k\X,T-mk) = - cog{(Pa - Pb) ■ X - {Ea - E^) ^ 

J V 2v^ 


X exp{ - {Ea - Eh) 


^(a^+nV/) (v.-Vb)2X2 


2 t ;2 


8((T^ -|- v‘^at)v‘^ 


} ( 21 ) 


where v = {v^ -fn/)/2. In terms of the average momentum P and energy E and their 

respective differences APand AE the interference term (21) ,to the leading order in {AP/P) 
and {AE/E), reduces to the form, 
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( 22 ) 


2 JC 

- cos{—(EAE — P ■ AP)} exp{—A} 


where the exponential damping factor, 


-( ) ^(;^) + (-^) ^ 


(23) 


and a2 = {a^ + (P/Efa^). 

In general neither AE nor AP is zero and they depend upon how the state is prepared. 
For example if the mommentum of the initial and the associated hnal state is measured 
then AP = 0. In the last section we have shown that for long distance propagation (s > 
m~^) P and E are not independent and are related by the mass shell constraint. The 
particular combination that appears in the phase difference turns out to be independent of 
the preparation and is hxed by the condition that each of the mass eigenstates be on shell, 


Ei^ — Pi^ = mj2 i = a,b 
=> 2EAE - 2P ■ AP = Am2 


(24) 


Substituting (24 ) in (22) we see that the interference term of (21) is given by , 

- cosi^X) e-^ , (25) 

The probability for the flavor conversion as a function of distance (20) is therefore, 
P{a^(3-X) = Y. - \Upa? \Uaa? 

a 


+ Y.- Wsa KJJk U„t,\cos(^XL^x + S) e-* (26) 

a^b 

where <5 = arg{Upa U*^aUpb Eab)- 

We see that the standard oscillation formula which was obtained for relativistic particles 
[1] is actually valid at all energies. In other words although the formula (26) is usually 
derived by taking the leading term in a {w?/P"^) series, our covariant calculation shows that 
there are actually no 0 {w?/P"^) corrections to (26). 

In the non-relativistic regime where P = {rria +mb)v/2, the standard kaon oscillation 
formula cosirUa — mb)T is recovered. 

There is no extra factor of two in the relativistic kaon oscillation formula as would have 
appeared without the wavepacket averaging . We see that although the two mass eigenstates 
have different proper times the wave packet overlap results in the average proper time 
appearing in the interference term. Therefore the interfernce phase is actually {ma — mb)s = 
{rua — mb){jna + mb)T/{Ea + Eb) = {nia — ml)T/2E. The average time prescription of 
[5] can therefore be justihed using the covariant propagator method. Other methods of 
showing that an extra factor of two does not appear in the flavor oscillations formula have 
been discussed in [6]. 
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Conversion without oscillation: The suppression factor A which goes with the oscillation 
term has some interesting new implications. If the parameters of the experiment are such 
that A becomes large then no oscillations is spacetime can be observed. What can be 
observed is a constant conversion probability. 

It has been noted ealier [8] one condition A must be small which implies that X must 
be smaller than the coherence length Lcoh = {2aP/AP) . We shall show below that the 
constraint Lose < Lcoh is not a sufficient condition to ensure the occurrence of flavor ocilla- 
tions in space. In the analysis of refs [8] only the uncertainty in the initial position was 
considered, in our analysis we have inclnded the contribution of at the nncertainty in time 
at which the neutrino was prodneed , which as it turns ont, makes the larger contribntion 
to the snppression term A. This happens when the nentrinos are prodneed from long lived 
resonances. The nentrino wavepacket which is prodneed has a spread in time with width at 
which cannot be smaller than the lifetime of the resonance whose decay prodnees it. The 
dominant contribntion in that case arises from the first term of A, 


A ~ {AEY ^ — 


AnrC r^2 

2V2E^ 


(27) 


where we have eqnated the initial time nncertainity at with r - the lifetime of the resonance 
that prodnees the nentrino, and we have ignored the spatial spread of the wave-packet a^ 
which in most experiments is many orders of magnitude smaller than at- In terms of the 
oscillation length Lose = ^tiE/{A nn?) the snppression factor A = (v^ttt/Losc)^- In order to 
observe oscillations the sonrce-detector distance L must be larger than Lose- Which that the 
minimnm source-detector distance in order that nentrino oscillations be observed is given 
by condition, Lmin = (7r)v^r. We shall show below that most nentrino experiments do 
not satisfy this criterion for observability of space oscillations. On the other hand spatial 
oscillations of flavor can be observed in the <h and B factories owing to the large width of 0 
and T resonances. 

For relativistic particles prodneed from long lived resonances the formula for the conver¬ 
sion probability which should be htted with the experiments is , 


P(z/q, —z/^) = - sin'^29 


,2.53Am‘^L^ 


cos( 


E 


exp{- 


A.79An?T 
■ E 


)^} 


( 28 ) 


where Am^ is the mass sqnare difference in ,L is the detector distance in m {km), r 
is the lifetime of the parent particle in the lab frame in m {km) and E is the energy in 
MeV {GeV). The limits on the valnes of Am^ and sin‘^26 obtained by htting the results of 
different experiments with the oscillation formula (28) are listed in Table I. 



TABLES 


TABLE 1. The asymptotic limits on Am? and sin?29 from different experiments according to 
the oscillation formula (28). r is the lifetime of the neutrino source in the lab frame, < > is the 

average u energy, L is the detector distance, Lmin is the minimum detector distance for observing 
oscillations in space and P is the experimental value of the conversion probability. 


Experiment (Source) 

r (m) 

< E^> (MeV) 

L (m) 

Lmin (m) 

1_L_ 1 

1 Arn^ 

sin^26 

LSND {^l) [9] 

658.6 

30 

30 

2927 


lyjBJitOliUi 

0.003 - 0.009 

LSND (tt) [18] 

17 

130 

30 

76 

|jj£gE3!B£IEIliMi 

1 

00 

o 

0.002 - 0.0082 

Karmen (tt) [11] 

7.8 

29.8 

17.5 

34.6 


< 0.16 


E776 (vr) [10] 

578 





< 0.2 


CCFR( K) [12] 






< 1.2 


Bugey(U,Pu) [13] 


1 5 

1 95 


1 < 0.75 X 10-" 

<Tcr^’ 

1 <0.15 
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FIG. 1. Lsnd allowed regions with the standard (between the dotted lines) and the covariant os¬ 
cillation formnlas (between the continnons line). The conversion probability formnlae are averaged 
over a Gaussian distribution of neutrino energy with mean SOMeV and width WMeV. 




Log(delta m^2 eV^2) 








-2.5 -2 -1.5 -1 -0.5 C 

Log(sin^2 2 theta) 


FIG. 3. The combined fit of the covariant oscillation formula with all experiments.The region 
between the continuous lines is allowed by the LSND fi experiment [9] .The region between the 
dotted lines is allowed by the LSND vr experiment [?] . Region ruled out by E776 [10] is above the 
dashed-dotted curve and by Karmen [11] is above dashed curve. The region above the top-most 
dashed curved is ruled out by CCFR [12]. 

12 




In the standard oscillation formula the oscillatory term averages to zero when Am^ > 
{AttE/L) . Using the wavepacket formula (28) however we see that the oscillatory term is 
is exponentially damped at much lower values of (provided r > (v^2L/47r)). If neu¬ 
trino sources have a large decay time which is the case in most experiments , the conversion 
probabilty is sensitive to mass differences Am^ > {\/22E/t). For this reason the LSND ex¬ 
periment which uses neutrinos from muon decay has two orders of magnitude more sensitive 
than KARMEN, BNL-E776 etc where the neutrinos are from tt and K decay. This is evident 
by comparing the £t of LSND fi DAR allowed region Figl with the Karmen allowed region 
Fig 2. In experiments where the source is in flight w.r.t the detector the gain in the source 
lifetime by the Lorentz factor 7 = {Et^/ttIt^) is somewhat compensated by the larger E^y in 
the denominator of A (27), so there is no gain in sensitivity with in flight neutrino sources 
(as in CCFR and BNL-E776 experiments). The reactor experiments where the neutrinos 
arise from the /5-decay of and nuclei , r ~ lO^s and therefore these experiments 

are in fact sensitive to values of Am^ as low as 10~^eV^. But the probability measurement 
in reactor experiments is poor compared to the accelerator experiments which is why they 
rule out only a small region of parameter space. If the recent LSND tt decay in flight result 
[18] is taken as an actual observation of > z/g conversion then it rules out most of the 
low Am^ region allowed by the LSND p decay experiment. 

Conclusions 

The covariant formulation gives a very different result for the oscillation probability 
formula compared to the standard treatment when the neutrino source has a lifetime larger 
than the baseline of the experiment. We use the condition that the minimum spread of the 
neutrino wave function along the time axis is given by the lifetime of the particle whose decay 
produces the neutrino. When neutrinos are produced from long lived particles like muons , 
the interference term in the conversion probability formula vanishes for much smaller values 
of the neutrino mass square difference compared with the standard formula. The reason this 
happens can be understood with the following picture. When the spread of the neutrino is 
large along the time axis, , its spread along the energy axis is small. The neutrino wave- 
functions for different masses are then energy eigenstates which are mutually orthogonal. 
The interference term which is the overlap of the wave-functions of different mass states 
therefore vanishes when these states become orthogonal. 

Acknowledgments I thank Terry Goldman, Walter Grimus and Marek Nowakowski for 
useful correspondence. 
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